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We report experimental results for the influence of a tilt angle (3 relative to gravity on 
turbulent Rayleigh-Benard convection of cylindrical samples. The measurements were 
made at Rayleigh numbers R up to 10 11 with two samples of height L equal to the 
diameter D (aspect ratio T = D/L ~ 1), one with L ~ 0.5 m (the "large" sample) and 
the other with L ~ 0.25 m (the "medium" sample). The fluid was water with a Prandtl 
number a = 4.38. 



In contrast to the experiences reported by Chilla et aL(2004) for a similar sample but 



with r ~ 0.5 {D = 0.5 and L = 1.0 m), we found no long relaxation times. 

For R — 9.4 x 10 10 we measured the Nusselt number TV as a function of [3 and obtained 
a small /3 dependence given by N{(3) = A/"o[l — (3.1 ± 0.1) x 10~ 2 |/3|] when (3 is in 
radian. This depression of N is about a factor of 50 smaller than the result found by 



Chilla et al.(2004)| for their V ~ 0.5 sample. 

We measured side- wall temperatures at eight equally spaced azimuthal locations on the 
horizontal mid-plane of the sample and used them to obtain cross-correlation functions 
between opposite azimuthal locations. The correlation functions had Gaussian peaks 
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centered about t\ c > that corresponded to half a turn-over time of the large-scale 

circulation (LSC) and yielded Reynolds numbers R^ c of the LSC. For the large sample 
and R = 9.4 x 10 10 we found Rf{0) = Rf{0) x [l + (1.85±0.21)|/3|-(5.9±1.7)/3 2 ]. Similar 
results were obtained also from the auto-correlation functions of individual thermometers. 
These results are consistent with measurements of the amplitude 6 of the azimuthal side- 
wall temperature- variation at the mid-plane that gave 8((3) = 6(0) x [1 + (1.84 ±0.45) |/3| — 
(3.1 ± 3.9)/3 2 ] for the same R. An important conclusion is that the increase of the speed 
(i.e. of R e ) with j3 of the LSC does not significantly influence the heat transport. Thus 
the heat transport must be determined primarily by the instability mechanism operative 
in the boundary layers, rather than by the rate at which "plumes" are carried away by 
the LSC. This mechanism apparently is independent of (3. 

Over the range 10 9 < R ^ 10 11 the enhancement of R% c at constant (3 due to the tilt 
could be described by a power law of R with an exponent of —1/6, consistent with a 
simple model that balances the additional buoyancy due to the tilt angle by the shear 
stress across the boundary layers. 

Even a small tilt angle dramatically suppressed the azimuthal meandering and the 
sudden reorientations characteristic of the LSC in a sample with (3 = 0. 

For large R the azimuthal mean of the temperature at the horizontal mid-plane differed 
significantly form the average of the top- and bottom-plate temperatures due to non- 
Boussinesq effects, but within our resolution was independent of (3. 



1. Introduction 

Turbulent convection in a fluid heated from below, known as Rayleigh-Benard convec- 
tion (RBC), has been under intense study for some time [for reviews, see e.g. Siggia(1994) 
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Kadanoff(2001)| Ahlers, Grossmann &: Lohse (2002)] . A central prediction of models for 



this system Kraichnan( 1 962 )[ Castaing et a/.(1998) Shraiman and Siggia (1990) Grossmann & Lohse (2001) 



is the heat transported by the fluid. It is usually described in terms of the Nusselt number 

_QL_ 

AXAT K ' 

where Q is the heat current, L the cell height, A the cross-sectional area, A the thermal 
conductivity, and AT the applied temperature difference. The Nusselt number depends 
on the Rayleigh number 

R = agATL 3 / ki> (1.2) 

and on the Prandtl number 

a = v I k . (1-3) 

Here a is the isobaric thermal expansion coefficient, g the acceleration of gravity, k the 
thermal diffusivity, and v the kinematic viscosity. 

An important feature of turbulent RBC is the existence of a large-scale circulation 



(LSC) of the fluid Krishnamurty & Howard(1981) . For cylindrical samples of aspect 



ratio r = L/D ~ 1 the LSC is known to consist of a single cell, with fluid rising 
along the wall at some azimuthal location 6 and descending along the wall at a lo- 
cation 9 + 7r [see, for instance, Qiu and Tong (2001a)| . As T decreases, the nature of 



the LSC is believed to change. For T < 0.5 it is expected Verzicco &: Camussi(2003)| 



Stringano & Verzicco(2005) Sun et al. (2005)| that the LSC consists of two or more con- 



vection cells, situated vertically one above the other. Regardless of the LSC structure, 
the heat transport in turbulent RBC is mediated by the emission of hot (cold) volumes 
of fluid known as "plumes" from a more or less quiescent boundary layer above (below) 
the bottom (top) plate. These plumes are swept away laterally by the LSC and rise (fall) 
primarily near the side wall. Their buoyancy helps to sustain the LSC. 



Guenter Ahlers, Eric Brown, and Alexei Nikolaenko 



In a recent paper Chilla et aL(2004) reported measurements using a cylindrical sample 
of water with a ~ 2.33 and with L = 1 m and D = 0.5 m for R ~ 10 12 . Their sample 
thus had an aspect ratio L ~ 0.5 at the borderline between a single-cell and a multi- 
cell LSC. They found exceptionally long relaxation times of JV that they attributed 
to a switching of the LSC structure between two states. Multi-stability was observed 



also in Nusselt-number measurements by Roche et al. (2004) for a T — 0.5 sample (see 



also Nikolaenko et al. (2005) for a discussion of these data). Chilla et al. also found 
that M was reduced by tilting the sample through an angle (3 relative to gravity by an 
amount given approximately by A/"(/?)/7V(0) ~ 1 — 2/3 when (3 is measured in radian. A 
reduction by two to five percent of Af (depending on R) due to a tilt by (3 ~ 0.035 of 



a T = 0.5 sample was reported as well recently by Sun et al. (2005) although in that 
paper the /3-dependence of this effect was not reported. Chilla et al. developed a simple 
model that yielded a depression of M for the two-cell structure that was consistent 
in size with their measurements. Their model also assumes that no depression of M 
should be found for a sample of aspect ratio near unity where the LSC is believed to 
consist of a single convection cell; they found some evidence to support this in the work 



of Belmonte et al. (1995) Indeed, recent measurements by Nikolaenko et al. (2005) for 
T = 1 gave the same N within 0.1 percent for a level sample and a sample tilted by 0.035 
rad. 

In this paper we report on a long-term study of RBC in a cylindrical sample with 
r ~ 1. As expected, we found no long relaxation times because the LSC is uniquely 
defined. The establishment of a statistically stationary state after a large change of R 
occurred remarkably quickly, within a couple of hours, and thereafter there were no 
further long-term drifts over periods of many days. 

We also studied the orientation 6q of the circulation plane of the LSC by measur- 
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ing the side- wall temperature at eight azimuthal locations Brown et al. (2005b)| . With 

the sample carefully leveled (i.e. (3 = 0) we found 6q to change erratically, with large 

fluctuations. There were occasional relatively rapid reorientations, as observed before 



by Sreenivasan et al. (2002) The reorientations usually consisted of relatively rapid ro- 
tations, and rarely were reversals involving the cessation of the LSC followed by its 
re-establishment with a new orientation. This LSC dynamics yielded a broad probability 
distribution-function P(9q), although a preferred orientation prevailed. When the sample 
was tilted relative to gravity through an angle /?, a well defined new orientation of the 
LSC circulation plane was established, P(9q) became much more narrow, and virtually 
all meandering and reorientation of the LSC was suppressed. 

We found that N was reduced very slightly by tilting the sample. We obtained Af((3) = 
A/o[l — (3.1 ± 0.1) x 10 2 1 /5 1 ] . This effect is about a factor of 50 smaller than the one 
observed by Chilla et al. for their T = 0.5 sample. 

From side-wall-temperature measurements at two opposite locations we determined 
time cross-correlation functions Ci_j. The C^j had a peak that could be fitted well by 
a Gaussian function, centered about a characteristic time t1 c that we interpreted as 
corresponding to the transit time needed by long-lived thermal disturbances to travel 
with the LSC from one side of the sample to the other, i.e. to half a turnover time of 
the LSC. We found that the /3-dependence of the corresponding Reynolds number R^ c 
is given by Rf{[3) = R c e c (0) x [1 + (1.85 ± 0.21) \(3\ - (5.9 ± 1.7)/3 2 ]. A similar result 
was obtained from the auto-correlation functions of individual thermometers. Thus there 
is an 0(1) effect of [3 on R e , and yet the effect of [3 on N was seen to be nearly two 
orders of magnitude smaller. We also determined the temperature amplitude 5 of the 
azimuthal temperature variation at the mid-plane. We expect 5 to be a monotonically 
increasing function of the speed of the LSC passing the mid-plane, i.e. of the Reynolds 
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number. We found 8{[3) = 6(0) x [1 + (1.84 ± 0.45) \/3\ - (3.1 ± 3.9)/3 2 ]. Thus, for small [3 

its /3-dependence is very similar to that of the Reynolds number. 

From the large effect of (3 on R e and the very small effect on Af we come to the im- 
portant conclusion that the heat transport in this system is not influenced significantly 
by the strength of the LSC. This heat transport thus must be determined primarily by 
the efficiency of instability mechanisms in the boundary layers. It seems reasonable that 
these mechanisms should be nearly independent of (3 when (3 is small. This result is con- 



sistent with prior measurements by Ciliberto et aL(1997)[ who studied the LSC and the 



Nusselt number in a sample with a rectangular cross section. They inserted vertical grids 
above (below) the bottom (top) plate that suppressed the LSC, and found that within 
their resolution of a percent or so the heat transport was unaltered. Their shadowgraph 
visualizations beautifully illustrate that the plumes are swept along laterally by the LSC 
when there are no grids and rise or fall vertically due to their buoyancy in the presence of 



the grids. Ciliberto et o/.(1997) also studied the effect of tilting their rectangular sample 
by an angle of 0.17 rad. Consistent with the very small effect of tilting on J\f found by 
us, they found that within their resolution the heat transport remained unaltered. 

We observed that the sudden reorientations of the LSC that are characteristic of the 
level sample are strongly suppressed by even a small tilt angle. 

2. Apparatus and Data Analysis 

For the present work we used the "large" and the "medium" sample and apparatus 



described in detail by Brown et al. (2005a) Copper top and bottom plates each contained 
five thermistors close to the copper-fluid interface. The bottom plate had imbedded 
in it a resistive heater capable of delivering up to 1.5 kW uniformly distributed over 
the plate. The top plate was cooled via temperature-controlled water circulating in a 
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Figure 1. Time evolution of the top- and bottom-plate temperatures and of the heat cur- 
rent for the large sample. Initially the temperature difference was AT ~ 0. At t = 0.6 hours 
the top- and bottom- plate regulators were given new set points corresponding to AT = 20° C 
(R — 9.43 x 10 10 ). (a): Top- (lower data set, T t ) and bottom- (upper data set, T b ) plate tem- 
peratures, (b): Heat current Q delivered by the temperature regulator designed to hold T& at a 
specified value, (c): Heat current Q on an expanded scale. The solid line represents a fit of an 
exponential function to the data for t > 1.2 h that gave a relaxation time tq = 0.48 h. 

double-spiral channel. For the Nusselt-number measurements a temperature set-point 
for a digital feedback regulator was specified. The regulator read one of the bottom-plate 
thermometers at time intervals of a few seconds and provided appropriate power to the 
heater. The top-plate temperature was determined by the temperature-controlled cooling 
water from two Neslab RTE740 refrigerated circulators. 

Each apparatus was mounted on a base plate that in turn was supported by three 
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legs consisting of long threaded rods passing vertically through the plate. The entire 

apparatus thus could be tilted by an angle (5 relative to the gravitational acceleration by 

turning one of the rods. The maximum tilt angle attainable was 0.12 (0.21) rad for the 

large (medium) sample. 

The Nusselt number was calculated using the temperatures recorded in each plate and 
the power dissipated in the bottom-plate heater. The side wall was plexiglas of thickness 
0.64 cm (0.32 cm) for the large (medium) sample. It determined the length L of the 
sample. Around a circumference the height was uniformly 50.62 ± 0.01 cm ( 24.76 ± 0.01 
cm) for the large (medium) sample. The inside diameter was D — 49.70 ± 0.01 cm ( 
D = 24.84 ± 0.01 cm) for the large (medium) sample. The end plates had anvils that 
protruded into the side wall, thus guaranteeing a circular cross section near the ends. 
For the large sample we made measurements of the outside diameter near the half-height 
after many months of measurements and found that this diameter varied around the 
circumference by less than 0.1%. 

Imbedded in the side wall and within 0.06 cm of the fluid-plexiglas interface were eight 
thermistors, equally spaced azimuthally and positioned vertically at half height of the 
sample. They yielded a relatively high (low) temperature reading at the angular positions 
where there was up-flow (down- flow) of the LSC. A fit of 

Ti = T c + 8 cos(in/4: - 9 ), i = 0,...,7 (2.1) 

yielded the mean center temperature T c , the angular orientation 9q of the LSC (relative 
to the location of thermistor 0), and a measure S of the LSC strength. 

We expect the size of S to be determined by the heat transport across a viscous 
boundary layer separating the LSC from the side wall. Thus 5 should be a monotonically 
increasing function of the LSC Reynolds number R e because the boundary-layer thick- 

1/2 

ness is expected to decrease with R e as l/R e , and because the azimuthal temperature 
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variation carried by the LSC near the boundary layer increases with R and thus with 

R e . However, the precise relationship between S and R e is not obvious. Experimentally 

we find, over the range 5 x 10 9 < R < 10 11 and for the large sample, that 8 is related 

to R by an effective power law 8 oc R 81 , whereas R e oc R 50 in this range, yielding 

8 oc i?g' 62 . We would then expect that 8 and R e will have a similar dependence on (3 (at 

least for small /?), albeit possibly with somewhat different coefficients. 

From time series of the Ti(t) taken at intervals of a few seconds and covering at least 

one day we determined the cross-correlation functions C l ,J (r) corresponding to signals 

at azimuthal positions displaced around the circle by 7r (i.e. j = i + 4). These functions 

are given by 

C^(t) = {[%{£) m))t] x m + r) - {T 3 (t)) t ]) t . (2.2) 

We also calculated the auto-correlation functions corresponding to i = j in Eq. 12.21 for 
all eight thermometers. 

Initially each sample was carefully leveled so that the tilt angle relative to gravity was 
less than 10 -3 radian. Later it was tilted deliberately to study the influence of a non-zero 
p on the heat transport. 

The fluid was water at 40°C where a = 3.88 x 1(T 4 K _1 , k = 1.53 x 1CT 3 cm 2 /s, and 
v = 6.69 x 1(T 3 cm 2 /s, yielding a = 4.38. 

3. The Nusselt number of a vertical sample 

3.1. Initial transients 

In Fig.^, we show the initial evolution of the top and bottom temperatures of the large 
sample in a typical experiment. Initially the heat current was near zero and Tf, and T t 
were close to 40°C. The sample had been equilibrated under these conditions for over 
one day. Near t = 0.6 h a new temperature set point of 50°C was specified for the bottom 
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Figure 2. The reduced Nusselt number N /R 1 ^ 3 for R = 9.43 x 10 10 for the large sample as 
a function of time during a single experimental run that lasted nine days. Each data point is 
based on temperature and heat-current averages over a time interval of two hours. The dotted 
line corresponds to the value estimated by using all the data. Note that the entire vertical axis 
covers only a change of 0.8%. The mean value is 0.06035, and the standard deviation from the 
mean is 5.1 x 10" 5 or 0.084%. 

plate, and the circulator for the top plate was set to provide T t ~ 30°C. From Fig. 
one sees that there were transients that lasted until about 0.9 h (1.2 h) for Tb (T t ). These 
transients are determined by the response time and power capability of the bottom-plate 
heater and the top-plate cooling water and are unrelated to hydrodynamic phenomena 
in the liquid. Figures ^ an d c show the evolution of the heat current. After the initial 
rapid rise until t ~ 0.8 h the current slowly evolved further to a statistically stationary 
value until t ~ 3 h. A fit of the exponential function Q(t) = Q x — AQ exp{~t / tq) to 
the data for t > 1.2 h is shown by the solid line in Fig.^ an d yielded a relaxation time 
tq = 0.48 ± 0.04 h. We attribute this transient to the evolution of the fluid flow. It is 
interesting to compare tq with intrinsic time scales of the system. The vertical thermal 
diffusion time t v = L 2 /k is 467 hours. Obviously it does not control the establishment 
of the stationary state. If we consider that it may be reduced by a factor of 1/Af with 
M = 263, we still obtain a time sale of 1.78 hours that is longer than tq. We believe 
that the relatively rapid equilibration is associated with the establishment of the top and 
bottom boundary layers that involve much shorter lengths l t and /{,. It also is necessary 
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Figure 3. The reduced Nusselt number M /R 1 ^ 3 as a function of the Rayleigh number R. 
Stars: data from |Chilla et al.(2004j| for a = 2.3 and T = 0.5. Solid circles: Data from 
|Nikolaenko et al. (2005")l for a = 4.38 and Y = 0.43. Open symbols: form |Nikolaenko et al. (2005) | 
for r = 0.67. Open squares: a — 5.42. Open circles: a = 4.38. Open diamonds: a — 3.62. 

for the large-scale circulation to establish itself; but, as we shall see, its precise Reynolds 
number is unimportant for the heat transport. In addition, the LSC can be created 
relatively fast since this is not a diffusive process. 

3.2. Results under statistically stationary conditions 

Figure ^ shows the behavior of the system only during the first six hours and does not 
exclude the slow transients reported by Chilla et al. that occurred over time periods of 
O(10 2 ) hours. Thus we show in Fig. [21 results for Af/R 1 / 3 from a run using the large 
sample that was continued under constant externally imposed conditions for nine days. 
Each point corresponds to a value of J\f based on a time average over two hours of the plate 
temperatures and the heat current. Note that the vertical range of the entire graph is 
only 0.8 %. Thus, within a small fraction of 1 %, the results are time independent. Indeed, 
during nearly a year of data acquisition for a T = 1 sample at various Rayleigh numbers, 
involving individual runs lasting from one to many days, we have never experienced long- 
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term drifts or changes of Af after the first few hours. This differs dramatically from the 

observations of Chilla et al. who found changes by about 2 % over about 4 days. We 

conclude that the slow transients observed by them for their T = 0.5 sample do not 

occur for r ~ 1. 

To document further the stationary nature of the system, we compared results from 
the large sample for Af obtained from many runs, each of one to ten days' duration, over 
a period of about five months [Nikolaenko et al. (2005) Funfschilling et al. (2005)] . The 
scatter of the data at a given R is only about 0.1%. This excellent reproducibility would 
not be expected if there were slow transients due to transitions between different states 
of the LSC. 

Although work in our laboratory with other aspect ratios has been less extensive, we 
also have not seen any evidence of drifts or transients for the larger V = 1.5,2,3, and 
6 [Funfschilling et al. (2005)) |Brown et al. (2005a)| nor for the smaller V = 0.67, 0.43, 
and 0.28 |Nikolaenko et al. (2005)1 |Brown et al - (2005a)| . It may be that T = 0.5, being 
near the borderline between a single-cell LSC and more complicated LSC structures 
Verzicco fc Camussi( 2003)1 |Stringano fc Verzicco(2005)| |Sun et al. (2005)] , s unique m 



this respect. 

In Fig.[21we compare results for Af/ i? 1 / 3 from our large sample [Nikolaenko et a/"72 005) 



with those reported by Chilla et al. (stars). Our results are larger by about 15%. To find 
a reason for this difference, we first look at the T and a dependence. The open (solid) cir- 
cles represent our data for a — 4.38 and T — 0.67 (0.43) and show that the dependence of 
Af on F is not very strong. The open squares (diamonds) are our results for F = 0.67 and 
cr = 5.42 (3.62) and indicate that Af actually increases slightly with a. Thus the lower 
values of Af (compared to ours) obtained by Chilla et al. for a = 2.3 and T = 0.5 can not 
be explained in terms of the T and a dependence of Af. Some of the difference can be 



Absence of slow transients, and the effect of imperfect vertical alignment, in turbulent Rayleigh-Benard convecti&fi 

1.004 




ca. 0.996 



0.99 



40 80 120 

time ( h ) 



160 



Figure 4. The reduced Nusselt number Af(/3)/Af(0) as a function of time for the large sample 
and R = 9.43 x 10 10 . Open circles: tilt angle [3 = 0. Solid circles: [3 = 0.087 radians. Open 
squares: [3 = 0.122 radians. All data are normalized by the average N[0) of the data for [3 = 0. 
Each data point is based on temperature and heat-current measurements over a two-hour period. 
The vertical dotted lines indicate the times when j3 was changed. 

attributed to non-Boussinesq effects that tend to reduce Af Funfschilling et al. (2005) . 
However, for the largest AT used by Chilla et al. (31°C) we expect this effect to be 
somewhat less than 1 % Funfschilling et al. (2005)] . Finally, the effect of the finite con- 
ductivity of the top and bottom plates comes to mind. This can reduce Af by several % 
when AT is large |Chaumat et a/.(2002)||Verzicco(2004)l|Brown et al. (2005a)| , but it is 
difficult to say precisely by how much. It seems unlikely that this effect can explain the 
entire difference, particularly at the smaller R (and thus AT) where it is relatively small. 



4. Tilt-angle dependence of the Nusselt number 

In Fig. 0] we show results for Af from the large sample at R — 9.43 x 10 10 . Each data 
point was obtained from a two-hour average of measurements of the various temperatures 
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and of Q. Three data sets, taken in temporal succession, for tilt angles = 0.000, 0.087, 

and 0.122 are shown. All data were normalized by the mean of the results for = 0. 

Typically, the standard deviation from the mean of the data at a given was 0.13%. 

The vertical dotted lines and the change in the data symbols show where was changed. 

One sees that tilting the cell caused a small but measurable reduction on N '. In Fig.[5]we 

show the mean value for each tilt angle, obtained from runs of at least a day's duration a 

each 0, as a function of \0\. One sees that Af decreases linearly with 0. A fit of a straight 

line to the data yielded 

Af{0) = Af [l - (3.1 ± 0.1) x 1O' 2 \0\] . (4.1) 

with Afo — 273.5. Simlar results for the medium cell are compared with the large-cell 
results in Fig. H3 At the smaller Rayleigh number of the medium sample the effect of 
on Af is somewhat less. Because the effect of on Af is so small, we did not make a more 
detailed investigation of its Rayleigh-number dependence. 

Chilla et al. proposed a model that predicts a significant tilt-angle effect on Af for 
T = 0.5 where they assume the existence of two LSC cells, one above the other. They 
also assumed that there would be no effect for T = 1 where there is only one LSC cell. 
Although we found an effect for our T = 1 sample, we note that it is a factor of about 
50 smaller than the effect observed by Chilla et al. for T = 0.5. 

5. Tilt-angle dependence of the large-scale circulation 

5.1. The orientation 

In Fig. |3 we show the angular orientation 9q (a) and the temperature amplitude S 
(b) of the LSC. For the first 8000 sec shown in the figure, the sample was level (0 = 
0.000 ±0.001). One sees that 9q varied irregularly with time. The probability-distribution 
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Figure 5. The Nusselt number M{f3) as a function of the tilt angle f3 for the large sample with 
R = 9.43 x 10 10 . Each point is the average over an entire run of duration one day or longer. 
Solid circles: (3 > 0. Open circles: j3 < 0. The solid line is a least-squares fit of a straight line 
(Eq.ETTJ to the data for (3 > 0. 

function P{Qq) is shown in Fig. [S] as solid dots. Essentially all angles are sampled by the 
flow, but there is a preferred direction close to 9q/2it = 0.6. At t = 8000 sec, the sample 
was tilted through an angle (3 — 0.087 radian. The direction of the tilt was chosen delib- 
erately so as to oppose the previously prevailing preferred orientation. As a consequence 
one sees a sharp transition with a change of 9q by approximately %. The temperature 
amplitude 5 on average increased slightly, and certainly remained non-zero. From this we 
conclude that the transition took place via rotation of the LSC, and not by cessation that 



would have involved a reduction of 5 to zero [see Brown et al. (2005b) . We note that 
9o(t) fluctuated much less after the tilt. The results for P(9q) after the tilt are shown in 
Fig. [5] as open circles. They confirm that the maximum was shifted close to 9q — 0, and 
that the distribution was much more narrow. 

In Fig. El we show P(9 ) for (3 = 0.122, (solid circles), 0.044 (open circles), and 0.026 
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Figure 6. The reduced Nusselt number M(P)/No as a fucntion of the tilt angle f3 for the large 
sample with R = 9.43 x 10 10 (solid circles) and the medium sample with R = 1.13 X 10 10 (open 
circles) . 

(solid squares). One sees that a reduction of (3 leads to a broadening of P(0q)- The square 
root of the variance of data like those in Fig. [5] is shown in Fig. ^| on a logarithmic scale 
as a function of [3 on a linear scale. Even a rather small tilt angle caused severe narrowing 
ofP(0 o ). 

5.2. The temperature amplitude 

In Fig.^Jwe show the temperature-amplitude 6((3) of the LSC as a function of (3. As was 
the case for A/", the data are averages over the duration of a run at a given (3 (typically 
a day or two). The solid (open) circles are for positive (negative) (3. The data can be 
represented well by either a linear or a quadratic equation. A least-squares fit yielded 



6(P) = 6(0) x [1 + (1.84 ± 0.45)|/3| - (3.1 ± 3.9)/3 2 



(5.1) 



with 6(0) = 0.164 K. 
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Figure 7. (a): The orientation 9o of the plane of circulation and (b): the temperature amplitude 
5 of the large-scale circulation as a function of time for the large sample and R — 9.43 x 10 10 . 
At t — 8000 sec the sample was tilted by an angle /3 = 0.087 radians relative to gravity. 

5.3. The Reynolds numbers 

Using Eg. 12. 21 we calculated the auto-correlation functions (AC) C % ' % , i = 0, . . . , 7, as well 
as the cross-correlation functions (CC) ,j = (i + 4)%8, i = 0, ... ,7, of the tempera- 
tures measured on opposite sides of the sample. Typical examples are shown in Fig. 
The CC has a characteristic peak that we associate with the passage of relatively hot or 
cold volumes of fluid at the thermometer locations. Such temperature cross-correlations 
have been shown, e.g. by Qiu and Tong (2001b) and Qiu and Tong (2002)} to yield delay 
times equal to those of velocity-correlation measurements, indicating that warm or cold 
fluid volumes travel with the LSC. The function 



C lJ (r) = -b exp 



7- 



b±exp 



(5.2) 



consisting of an exponentially decaying background (that we associate with the random 
time evolution of #o) and a Gaussian peak, was fitted to the data for the CC. The fitted 
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Figure 8. The probability distribution P(6o) of the orientation 9o of the plane of circulation of 
the large-scale flow for the large sample and R = 9.43 x 10 10 as a function of do- Solid circles: 
tilt angle (3 = ± 0.001. Open circles: (3 = 0.087 radians. 
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Figure 9. The probability distribution P(6o) of the orientation 9o of the plane of circulation of 
the large-scale flow for the large sample and R = 9.43 x 10 10 at three tilt angles (3. Solid circles: 
(3 = 0.122. Open circles: (3 = 0.044. Solid squares: (3 = 0.026. 

function is shown in Fig. 1121 as a solid line over the range of r used in the fit. It is an 
excellent representation of the data and yields the half turnover time T/2 = t\ J of the 
LSC. Similarly, we fitted the function 

r / v 2i r / . . s 

T — to \ 



C 1 ' 1 {t) = b exp 



biexp 











[-(?)' 


+ b 2 exp 







(5.3) 



Absence of slow transients, and the effect of imperfect vertical alignment, in turbulent Rayleigh-Benard convecti&u 







• 






> 










•• 


• 






• 


CD 


• 


• 




• • 





-0.05 0.00 0.05 0.10 



p ( rad ) 

Figure 10. The square root of the variance oe of the probability distribution P(6o) of the orien- 
tation 9 of the plane of circulation of the large-scale flow for the large sample and R — 9.43 x 10 10 
as a function of the tilt angle /3. 



0.20 




0.161 1 1 1 1 1 1— 

0.00 0.04 0.08 0.12 

Ipl(rad) 



Figure 11. The time-averaged temperature amplitude S((3) of the LSC for the large sample as 
a function of \f3\. Solid circles: (3 > 0. Open circles: (3 < 0. For this example R = 9.43 x 10 10 . 

to the AC data. It consists of two Gaussian peaks, one centered at t = and the other 
at r = t l 2 , and the exponential background. We interpret the location t£ of the second 
Gaussian peak to correspond to a complete turnover time T of the LSC. 
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FlGURE 12. The cross-correlation lunction C 04 (r) between thermometers and 4 and the 
auto-correlation function C 00 (r) of thermometer for the large sample. The solid lines are 
fits of Eas. 15.21 and 15.31 to the data. They also indicate the range of r used for the fits. For this 
example the tilt angle was = -0.009 and R = 9.43 x 10 10 . 

In terms of the averages < t{ > and < ft} 1 > over all 8 thermometers or thermometer- 
pair combinations we define Qiu and Tong (2002) Grossmann fc Lohse (200'2)] the Reynolds 
numbers 



{Lf <t\* >)[L/ V ) 



(5.4) 



and 



RT = (2L/ < >){L/v) . (5.5) 

Here the length scale 2L was used to convert the turnover time T into a LSC speed 
2L/T. For r = 1, the length AL might have been used instead, as was done for instance 



by Lam et al. (2002) This would have led to a Reynolds number larger by a factor of two. 
In Fig. 113b and b we show i?g C (|/3|) and R^ c (\/3\) respectively. The solid circles are for 
positive and the open ones for negative (3. One sees that i?g C (|/3|) and R^ c (\(3\) initially 
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Figure 13. (a): The Reynolds number R c e c (\f3\) obtained from the temperature cross-correlation 
functions, and (b): the Reynolds number -Re C (|/3|) obtained from the temperature 
auto-correlation functions, of the LSC for the large sample and R = 9.43 x 10 10 as a func- 
tion of the absolute value \/3\ of the tilt angle. Solid circles: /? > 0. Open circles: < 0. For this 
example R = 9.43 x 10 10 . 

grow linearly with 0, but the data also reveal some curvature as \0\ becomes larger. Thus 
we fitted quadratic equations to the data and obtained 



Rf{0) = Rf(Q) x [1 + (1.85 ± 0.21)|/3| - (5.9 ± l.7)0 2 } 



(5.6) 



and 



R a e c (0) = R a e c (Q) x [1 + (1.72 ± 0.38)|/3| - (4.1 ± 3.2)/3 2 ] (5.7) 
with Rf{Q) = 10467 ± 43 and R a e c {0) = 10565 ± 82 (all parameter errors are 67% conn- 
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Figure 14. The Reynolds-number ratio R c £ c (ft) / R c £ c (Q) (solid circles) obtained from the temper- 
ature cross-correlation functions, and the amplitude ratio 5(f3)/5(0) (open circles), of the LSC 
of the medium sample for R = 1.13 x 10 10 as a function of f3. 

dence limits). The results for i?g C (0) and i?g C (0) are about 10% higher than the prediction 



by Grossmann & Lohse (2002) for our a and R. The excellent agreement between Rf 
and _R" C is consistent with the idea that the CC yields T/2 and that the AC gives T. As 
expected (see Sect.EJ, the /3-dependences of both Reynolds numbers are the same within 
their uncertainties. It is interesting to see that the coefficients of the linear term also 
agree with the corresponding coefficient for 5 (Eq. 15.11) . This suggests that there may be 
a closer relationship between 6 and R e than we would have expected a priori. However, 
the coefficient of the linear term in Eq. 15.61 or 15.71 is larger by a factor of about 50 than 
the corresponding coefficient for the Nusselt number in Eq. 14.11 

In Fig. El we show measurements of R c e c and of 5, each normalized by its value at 
(3 = 0, as a function of (3 for the medium sample and R = 1.13 x 10 10 . For this sample 
we were able to attain larger values of (3 than for the large one. One sees that 6 and Rf 
have about the same j3 dependence for small (3, but that 5 then increases more rapidly 
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Figure 15. The Reynolds-number ratio Rf((3 = 0.122)/iC(0) of the LSC as a function of R. 
Open circles: medium sample. Solid circles: large sample. Solid line: powerlaw with an exponent 
of 1/6. 

than i?g C as (3 becomes large. Although we do not know the reason for this behavior, it 
suggests that the larger speed of the LSC enhances the thermal contact between the side 
wall and the fluid interior. 

The Rayleigh-number dependence of R c e c at constant (3 is shown in Fig. Here the 
open (solid) circles are from the medium (small) sample. There is consistency between 
the two samples, and the data can be described by a power law with a small negative 
exponent. The solid line is drawn to correspond to an exponent of —1/6. 

5.4. A model for the enhancement of the Reynolds number 

As seen in Fig. [H] the LSC assumes an orientation for which gravity enhances the velocity 
above (below) the bottom (top plate), i.e. the LSC flows "uphill" at the bottom where it 
is relatively warm and "downhill" at the top where it is relatively cold. This leads to an 
enhancement of the Reynolds number of the LSC. As suggested by Chilla et al. (2004)) 
one can model this effect by considering the buoyancy force per unit area parallel to 
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the plates. This force can be estimated to be plg/3aAT/2 where I is the boundary-layer 

thickness. It is opposed by the increase of the viscous shear stress across the boundary 

layer that may be represented by pvvl /I where u! is the extra speed gained by the LSC 

due to the tilt. Equating the two, substituting 

I = L/{2N) , (5.8) 

solving for vl , using Eg. 11.21 for R, and defining R' e = (L/v)u' one obtains 

for the enhancement of the Reynolds number of the LSC. From our measurements at 
large R we found that R e Ahlers et al. (2005)] and N Nikolaenko et al. (2005) can be 



represented within experimental uncertainty by 

R e = 0.0345i? 1/2 , (5.10) 

TV = 0.0602i? 1/3 , (5.11) 

giving 

^ = 1.00 x 10 3 R^ 1 / 6 a~ 1 p . (5.12) 

Re 

For our a = 4.38 and R = 9.43 x 10 10 one finds R' e /R e = 3.4/3, compared to the 
experimental value (1.9 ± 0.2)/? from Rf [Eq.[5T5| and (1.7 ± 0.4)/? from R^ c [Ea.K7\. 
We note that the coefficient 1.00 x 10 3 in Eq. 15.121 depends on the definition of R e 
given in Eqs. 15.41 an 15.51 that was used in deriving the result Eq. 15.101 If the length 
scale AL had been used instead of 2L to define the speed of the LSC, as was done for 
instance by Lam et a/.(2002)[ this coefficient would have been smaller by a factor of 



two, yielding near-perfect agreement with the measurements. In Fig. 1151 one sees that 
the predicted dependence on i? -1 / 6 also is in excellent agreement with the experimental 
results. However, such good agreement may be somewhat fortuitous, considering the 
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approximations that were made in the model. Particularly the use of Eq. 15.81 for the 

boundary-layer thickness is called into question at a quantitative level by measurements 



of Lui & Xia(1998) that revealed a significant variation of I with lateral position. In 



addition, it is not obvious that the thermal boundary-layer thickness I should be used, 



as suggested by Chilla et aL(2004) to estimate the shear stress; perhaps the thickness of 



the viscous BL would be more appropriate. 



In discussing their V = 0.5 sample, Chilla et ai. (2004) took the additional step of 



assuming that the relative change due to a finite (3 of Af is equal to the relative change 
of R e . For our sample with T = 1 this assumption does not hold. As we saw above, the 
relative change of Af is a factor of about 50 less than the relative change of R e . The 
origin of the (small) depression of the Nusselt number is not so obvious. Naively one 
might replace g in the definition of the Rayleigh number by g cos([3); but this would 
lead to a correction of order (3 2 whereas the experiment shows that the correction is of 
order f3, albeit with a coefficient that is smaller than of order one. The linear dependence 
suggests that the effect of (3 on Af may be provoked by the change of R e with [3, but not 
in a direct causal relationship. 

6. Tilt-angle dependence of reorientations of the large-scale 
circulation 

It is known from direct numerical simulation 



Hansen et a/.(1991)| and from several ex- 



periments Cioni et a/.(1997) Niemela et al. (2001)| Sreenivasan et al. (2002)||Brown et al. (2005b) 



that the LSC can undergo relatively sudden reorientations. Not unexpectedly, we find 
that the tilt angle strongly influences the frequency of such events. For a level sam- 
ple {(3 = 0) we demonstrated elsewhere Brown et al. (2005b) that reorientations can 



involve changes of the orientation of the plane of circulation of the LSC through any 
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Figure 16. The number of reorientation events per hour of the angular orientation of the 
plane of circulation of the LSC in the large sample for R — 9.43 x 10 10 . 

angular increment A9, with the probability P(A8) increasing with decreasing AO. Thus, 
in order to define a "reorientation" , we established certain criteria. We required that 
the magnitude of the net angular change \A8\ had to be greater than A0 min = (2ir)/8. 
In addition we specified that the magnitude of the net average azimuthal rotation rate 
\8\ = \A8/At\ had to be greater than 8 m in = 0.1/T where T is the LSC turnover time 



and At is the duration of the reorientation (we refer to Brown et al. (2005b) for further 
details). Using these criteria, we found that the number of reorientation events n(J3) at 
constant R = 9.43 x 10 10 decreased rapidly with increasing |/3|. These results are shown 
in Fig. It is worth noting that nearly all of these events are rotations of the LSC 
and very few involved a cessation of the circulation. A least-squares fit of the Gaussian 
function 

n{J3) = N exp[-{p - (3 ) 2 /w 2 ] (6.1) 
to the data yielded N = 1.23 ± 0.06 events per hour, /3 = 0.0093 ± 0.0010 rad, and 
w = 0.0251 ± 0.0015 rad. It is shown by the solid line in the figure. 

We note that the distribution function is not centered on (3 — 0. The displacement of 
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the center by about 9 mrad is much more than the probable error of (3. We believe that 

it is caused by the effect of the Coriolis force on the LSC that will be discussed in more 

detail elsewhere Brown et al. (2005c)| . 



7. Tilt-angle dependence of the center temperature 

We saw from Fig. ^2 that the increase of R e with (3 led to an increase of the amplitude 
8 of the azimuthal temperature variation at the horizontal mid-plane. An additional 
question is whether the tilt-angle effect on this system has an asymmetry between the top 
and bottom that would lead to a change of the mean center temperature T c (see Ea. l2.1|) . 



Chilla et a/.(2004) report such an effect for their T = 0.5 sample. For a Boussinesq sample 
with [3 = we expect that T c = T m with T m = (T t + T b )/2 (T f and T b are the top 
and bottom temperatures respectively), or equivalently that A t = T c — T t is equal to 
At = T b — T c . A difference between Af, and A t will occur when the fluid properties 
have a significant temperature dependence |Wu fc Libchaber(1991)[ Zhang et a/.(1997)] , 



i.e. when there are significant deviations from the Boussinesq approximation. For the 
sequence of measurements with the large apparatus and R = 9.43 x 10 10 as a function of 
[3 the mean value of AT = T b —T t was 19.808±0.018°C and T c -T m was 0.97°C, indicating 
a significant non-Boussinesq effect. In Fig. 117b , we show A t and Af, as a function of (3. 
One sees that increasing [3 does not have a significant effect for our T = 1 sample. This 
is shown with greater resolution in Fig. 117b where T c — T m — (A t — Af,)/2 is shown. 
We believe that the small variation, over a range of about 5 x 10~ 3o C, is within possible 
systematic experimental errors and consistent with the absence of a tilt-angle effect. 
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Figure 17. (a): The temperature difference in Kelvin between the bottom and the center 
(A;, = Tb — T c , soiid circles) and the center and the top (At = T c — T t , open circles) as a 
function of the tilt angle /3. (b) : The temperature difference T c — T m = (At — A(,)/2 between the 
center temperature T c and the mean temperature T m = (T t + Tt)/2. The center temperature is 
the average of the values given by the eight side-wall thermometers. 

8. Summary 

In this paper we reported on an experimental investigation of the influence on turbulent 
convection of a small tilt angle (3 relative to gravity of the axes of two cylindrical Rayleigh- 
Benard samples. The aspect ratios were F ~ 1. 

Where there was overlap, there were significant differences between our results and 



those obtained by Chilla et a/.(2004) for a T = 0.5 sample. We found our system to 
establish a statistically stationary state quickly within a couple of hours, after a Rayleigh- 



number change whereas Chilla et aL(2004) found long transients that they attributed to 
changes of the LSC structure. We found a very small depression of the Nusselt number 
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N with increasing (3, by about 4% per radian at small (3. Chilla et aL(2004) found a 
decrease by 200% per radian for their sample. 

In contrast to the very small effect of (3 on Af, we found an increase of the Reynolds 
number R e by about 180% per radian for small (3. The small effect on Af in the presence of 
this large change of R e indicates that the heat transport does not depend strongly on the 
speed of the LSC sweeping over the boundary layers. Instead, Af must be determined by 
instability mechanisms of the boundary layers, and the associated efficiency of the ejection 
of hot (cold) volumes (so-called "plumes") of fluid from the bottom (top) boundary layer. 

It is interesting to note that the strong dependence of R e on (3 in the presence of only 
a very weak dependence of Af on [3 can be accommodated quite well within the model 



of Grossmann & Lohse (2002) The Reynolds number can be changed by introducing a 
/3-dependence of the parameter a(f3) in their Eqs. (4) and (6). As pointed out by them, 
a change of a has no influence on the predicted value for Af. 

We also measured the frequency of rapid LSC reorientations that are known to occur 
for (3 = 0. We found that such events are strongly suppressed by a finite (3. Even a 
mild breaking or the rotational invariance, corresponding to j3 ~ 0.04, suppresses re- 
orientations almost completely. 
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